A superconducting qubit with Purcell protection and tunable coupling 
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We present a superconducting qubit for the circuit quantum electrodynamics architecture that has 
a tunable coupling strength g. We show that this coupling strength can be tuned from zero to values 
that are comparable with other superconducting qubits. At g = the qubit is in a decoherence 
free subspace with respect to spontaneous emission induced by the Purcell effect. Furthermore we 
show that in the decoherence free subspace the state of the qubit can still be measured by either a 
dispersive shift on the resonance frequency of the resonator or by a cycling- type measurement. 



Quantum decoherence is one of the major problems 
facing quantum information processing. To overcome 
this problem the theories of quantum error correction [1] 
and decoherence free subspaces (DFS) [2] have been de- 
veloped. A DFS is a subspace of a system which exploits 
symmetries in the decoherence process to allow the sys- 
tem to be completely decoupled from the environment. 
As an example, the spontaneous decay of a multilevel 
atom into the same bath can be cancelled for one of the 
states by quantum interference [3]. 

In recent years superconducting qubits have emerged 
as candidates for quantum information processing [4]. 
These are systems which are designed using Josephson 
junctions to make low loss non-linear oscillators. They 
are designed so that two levels (qubit) can be isolated, 
controlled and measured, properties which are usually 
mutually exclusive. With sweet spot operations [5-7] 
and material engineering [8] there has been tremendous 
progress. This is evidenced by the recent demonstration 
of two qubit quantum algorithms [9], high fidelity sin- 
gle qubit gates [10], high fidelity two [11, 12] and three 
[13, 14] qubit entangled states, and Bell violation [15]. 

Currently the most successful superconducting qubits 
are the flux [6] , phase [8] , and transmon [7] as these qubits 
are essentially immune to offset charge (charge noise) by 
design. The transmon receives its charge noise immunity 
by operating at a point in parameter space where the 
energy level variations with offset charge are exponen- 
tially suppressed. This suppression has experimentally 
been observed and resulted in these qubits being approx- 
imately T\ limited (T2 ~ 2Ti) in the circuit quantum 
electrodynamics (QED) architecture [16]. In this archi- 
tecture the qubits are coupled to a coplanar waveguide 
resonator through a Jaynes-Cummings Hamiltonian op- 
erated in the dispersive regime [17]. This resonator acts 
as the channel to control, couple, and readout the state 
of the qubit (see Fig. 1 A). 

In the circuit QED architecture, a significant source of 
T\ has been shown to be Purcell decay [18]. This is a fun- 
damental relaxation that arises when a qubit is coupled 
to a resonator. It can be understood as dressing of the 




FIG. 1. (color online) A) is a schematic of the circuit QED 
architecture. B) is a schematic of the proposed three island 
device. The islands are connected by SQUIDs and the arrows 
are used to indicated the dipole and quadrupole moments of 
the device. C) Is the circuit model of the device and the 
variables are explained within the text. 



qubit, the dressed qubit contains a photonic part propor- 
tional to the coupling strength g divided by the detuning 
A between the qubit and resonator. This photonic part 
of the qubit will relax at a rate given by the decay rate 
of the resonator n and as a result the dressed qubit will 
relax at a rate given jf ur = ng 2 / 'A 2 . To overcome the 
Purcell decay we can either work at large detunings, use 
a Purcell filter [19], or design a Purcell protected qubit. 

In this letter we will present a three island device that 
has the properties of a qubit (two levels, arbitrary con- 
trol and measurement) and has the ability to indepen- 
dently tune both the resonance frequency and coupling 
strength g, whilst still exhibiting exponential suppres- 
sion of the charge noise and maintaining an anharmonic- 
ity equivalent to that of the transmon and phase qubit. 
This tunable coupling qubit (TCQ) can be tuned from a 
configuration which is totally Purcell protected from the 
resonator g = (in a DFS) to a position which couples 
strongly to the resonator with values comparable to those 
realized for the transmon. Furthermore, we show that 
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in the DFS position a strong measurement can be per- 
formed. The TCQ only needs to be moved from the DFS 
position when single and two qubit gates are required and 
as such in the off position all multi-qubit coupling rates 
are zero. That is, the TCQ in the circuit QED archi- 
tecture (see Fig. 1 A) is its own tunable coupler to any 
other TCQ, going beyond the nearest neighbour tunable 
couplers presented in Refs. [20] and [21]. 

The essential idea behind the TCQ is a three island 
version of the transmon as shown in Fig. 1 B. This de- 
vice, like the transmon, has a dipole moment between 
each island. These dipoles can add in parallel, result- 
ing in a larger dipole moment, or in antiparallel, creating 
a quadrupole moment. This device now supports two 
different modes and with only the dipole moment be- 
ing able to couple to the resonator. That is, the mode 
corresponding to the quadrupole moment is a DFS with 
respect to Purcell decay and can be used for the storing 
of quantum information. Due to the capacitance between 
the top and lower island Cj these modes couple and the 
ratio of quadrupole to dipole moment can be tuned by 
changing the energy of the upper and lower island. 

The reduced circuit model we use for this device is 
shown in Fig. 1 C. The solid green lines represent the 
components associated with the TCQ with light indicat- 
ing the upper " + " and lower island " — " and the dark 
representing the center island. Each island is connected 
by a Josephson junction of energy Ej ± and capacitance 
C± which is taken to include the Josephson capacitance. 
The resonator is approximated by a parallel LC circuit 
(blue dashed lines) with inductance L r and capacitance 
C r . C c represents the capacitive interaction between the 
TCQ and resonator. Finally the dotted red lines repre- 
sent charge noise resulting from voltage fluctuations V g± 
that occur inside the device (C g± represent the capacitor 
coupling for these fluctuations). Note a similar system 
was presented in Ref. [23] for observation of giant non- 
linear Kerr effects in circuit QED. 

The Hamiltonian is obtained using the method out- 
lined in Ref. [22] with the three degrees of freedom be- 
ing the phase across the junctions 7 + = (0 + — 0_)/<I>o, 
7_ = 27r0_/ < l ) o, and the flux across the resonator <j) r . 
Here $0 = h/2eis the flux quantum and <p± are the node 
fluxes defined in Fig. 1 C. We find H = H T + Hi + H R 
where = huo r a^a is the Hamiltonian of the resonator 
with uj r = 1/ ^L r C' r and a being the standard annihi- 
lation operator, Ht the Hamiltonian for the TCQ given 
by 

H T = Y^^ E c ± (n±-n' g± ) 2 - y ^E J± cos( 7± )+4£; / n + n_ 
± ± 

(i) 

with charging energy Ec ± = e 2 /2C±, interaction en- 
ergy Ej = e 2 /C}, and dimensionless gate voltage n' g± = 
n g± + n g C^/Cj where n g± = C g± V g± /2e. Note the 
prime above the capacitance indicate that they have been 




0.0 0.2 0.4 0.6 0.8 1.0 20 40 60 80 



Ej/ E c Ej/ E c 

FIG. 2. (color online) A) Eigenenergies of the TCQ Hamilto- 
nian as a function of Ej/Ec for Ej ± = 50Ec- Solid lines are 
from a numerical diagonalization and dashed lines are from 
the coupled anharmonic oscillator model. B) Charge disper- 
sion |£ q ( m )| as a function of the ratio Ej/Ec for Ei = —Ec 
(solid lines) and Ei = (dashed lines). 



renormalized by the interactions. Finally Hi represents 
the interaction of the resonator with the TCQ and is 

Hi = 2e 2 V rms (/3+n+ + /3_n_) (-m f + ia) (2) 

where /3± = C c C^ t /[C^ + C^_ + (Cj + C C )(C E+ + C E J] 
with C s± = C± + C g± + C c , and V rms = y / huj r /2C' r . In 
the limit where the TCQ is isolated (C c = 0) and sym- 
metric (drop all ± dependence in capacitors) then Eq + = 
E c _ = E c with E c = e 2 (d + C E )/2(C| + 2C I C Yj ) and 
Ej = —2EcCj/{Cj + C$]). Ej can be tuned from zero 
to —2Ec by modifying Cj, governed by the direct ca- 
pacitance between the upper and lower island. This can 
be made much larger then the interaction energy between 
two superconducting qubits that are coupled virtually by 
a resonator [17]. The eigenenergies of the TCQ Hamil- 
tonian are shown in Fig. 2 A (solid lines) as a function 
of EjjE c for E J± = 50 E C - The system has a "V" like 
structure with two levels in the first excitation manifold 
and three in the next. As we increase Ej the degenera- 
cies in the manifolds are lifted and we have a multilevel 
atom with non ladder like structure. 

In the limit of large Ej ± /Ec± the values of j± are 
restricted to be around zero; thus we can expand the 
cosine terms up to fourth order [7]. That is, we can 
model the system as two coupled anharmonic oscillators 
with the Hamiltonian Hr eff = fl+ + H_ + H c with H± = 
h[uj± + 5± (b ] ± b±-l)/2]b ] ± b± and H c = HJ(b+bl_ + b\b_) 
with u± = yj SEj ± Ec ± I h — Ec±/h, S± = —Ec±/h, and 
J = Ei(Ej + Ej_/E c+ Ec_) 1/4 /V2h. This Hamiltonian 
can be approximately diagonalized (in orders of 5±/(oj+ — 
Cj-) with the transformation D = exp[A(6 + &L — &+&-)] 
with A = tan _1 (2J/^)/2 + and 77 = cj + — uo- — {5+ — 
5-)/2. Here = for (J, 77) > 0, = tt/2 for 77 < 0, and 
6 = 7r for (— J, 77) > 0. Doing this we find 

Ht m = ft J>± + S -^(blb± - l)P ± b ± + S c b\_b + blb_, 
± 1 

(3) 
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where uj± = ou± + (6± - S±)/2 + (£+ + 5_)J 2 /2fi 2 ± 
/i/2 T r//2, 5± = ((5 + + (5_)(l + r ? 2 //i 2 )/4±^ + -(5_)/2 / i 
and £ c = 2 J 2 (<J+ + £_)//i 2 with /i = ^4 J 2 + and the 
tilde indicating the diagonalized frame. The coupling has 
induced a conditional anharmonicity S c , it is this anhar- 
monicity that makes this system different to two coupled 
qubits, it ensures that En is not equal to Eqi + £ao- 
Here we have introduced the notation that superscript 
ij refers to i excitations in the dark mode ( " + " ) and 
j excitations in the bright mode (" — "). The choice of 
these names will become clearer latter. The dotted lines 
in Fig. 2 A are the predictions from this effective model, 
which agree well with the full numerics. Thus from the 
effective model the anharmonicities are all around Eq 
provided \J\ > \rj\. That is the TCQ has not lost any 
anharmonicity in comparison to the transmon or phase 
qubit and with simple pulse shaping techniques arbitrary 
control of the lowest three levels will be possible [24] . We 
will now introduce the notation that the qubit is formed 
by the space |0) = |00), |1) = |10) and \m) = |01) is the 
measurement state. 

Since charge fluctuations are one of the leading sources 
of noise in superconducting circuits we want to ensure 
that quantum information in the TCQ is not destroyed 
by charge noise. Following Ref. [7], the dephasing time 
for the qubit and m level will scale as l/|£ g ( m ) | where 
£ q (jn) is the the peak to peak value for the charge dis- 
persion of the — 1 and — m transition respectively. 
The dispersion in the energy levels arises from the gate 
charges n Qa and the fact the the potential is periodic. 
This can not be predicted with the coupled anharmonic 
oscillator model and as such is investigated numerically. 
We expect, that like the transmon, this will exponentially 
decrease with the ratio of Ej/Ec as in this limit the ef- 
fects of tunneling from one minima to the next becomes 
exponentially suppressed. This is confirmed in Fig. 2 
B where the have plotted \£ q (rn)\/ Eq(m) (the numerical 
maximum and minimum of the energy level over n' ga ) as 
a function of Ej/E c for Ej = E c and Ej = (trans- 
mon limit). That is, the TCQ has the same charge noise 
immunity as the transmon. 

We now investigate how the two modes of the TCQ 
couple to the resonator. We start by applying the anhar- 
monic oscillator model to Eq. (2). Doing this we find that 
= hT,±9±(ab t ± + a^b±) where g± = #± cos(A) ± 
# T sin(A) with g± = V2e 2 f3 ± V rms (E J± /8E c± ) 1/4 /h. 
Thus one of these coupling strengths can be set to zero 
for an appropriately chosen A. For the case when J is 
negative it is the g+ rates which can be set to zero, 
hence the name dark for the " + " and bright for the 
" — " state. We have numerically confirmed that this 
model approximately predicts the matrix elements (cou- 
pling rates) for the first two manifolds (six levels). Us- 
ing these six levels we can use perturbation theory to 
find an effective description of the situation where both 



the bright and dark modes of the TCQ have a disper- 
sive interaction with the resonator, |A±| ^> \g±\ with 
A± = uj± — uj r . To second order in g±/A± we find that 
the resonator-TCQ interaction induces both a Lamb shift 
on the TCQ and a TCQ state dependent pull on the res- 
onator. The qubit (0 — 1) transition frequency becomes 
uo q = Cj+ + g 2 _/A + and the — m transition frequency 
becomes uj m = cD_ + g 2 _ /A_. The resonator frequency 
is uj r + Xk with resonator pull xo = — g 2 -/ A_ — < g 2 _/A + , 
Xi = (5+ ~ A+)^ 2 /A + (£ + + A+) - g 2 _/(5 c + A_), and 
Xm = (S- ~ A_)£/A_(L + A_) - ~g\ /& + A+). All 
states |0), |1), and \m) have a different resonator fre- 
quency and thus can be measured by probing the trans- 
mission of the resonator. 

To calculate the induced decay rate on the qubit by 
relaxation of the resonator (Purcell decay) we can use 
the same perturbation theory. This amounts to evalu- 
ating 7^ ur = ft|(l|a|0)| 2 where the bar is the first order 
correction to the TCQ levels from the resonator-TCQ 
interaction. Doing this we find j p = Kg^/A 2 ^. Thus, 
setting g+ to zero the Purcell effect (to second order) 
will be canceled and yet the difference between the res- 
onator pull for the |0) and |1) state is non-zero and 
given by x = 5 c g 2 _/A-(A- + S c ). This has the same 
functional form as the transmon (when \J\ > \rj\) and 
is detectable with current microwave electronics. Fur- 
thermore, signal-to-noise (SNR) arguments from Ref. [25] 
carry over to this system and we find SNR = 4nf]xT m . 
Here rj is efficiency of collecting the photons emitted 
from the resonator, n is the number of the photons in 
the resonator which should not exceed n cr i t = A 2 _/4^ 2 _ 
[25], and T m is the measurement time. Taking realis- 
tic values T m = 1/is, x/2tt = 10 MHz, and rj = 1/20 
gives a SNR around 13n. However, much higher val- 
ues can be obtained by using the following protocol: 
Set uj m = uj r and have C c large enough to ensure that 
\uj q — uj r \ ^> \g±\ then the — m transition will vacuum 
Rabi split the cavity transmission if the qubit is in the 
state otherwise the transmission will be the bare res- 
onator. That is, the heterodyne power ft|(a)| 2 in steady 
state will be ^ 2 (g 2 _ + K 2 /l6)/(g 2 _ + ^ 2 /16 + £ 2 ) 2 [26] 
for the |0) and 4£ 2 /ft for the |1) state with £ being the 
cavity drive amplitude at drive frequency co r . Taking 
SNR = rjT m K>(\(a) o| 2 — |(&)i| 2 ) we find values as large as 
1000 for g- /2n = 100 MHz, n/2n = 10 MHz, and £ = <?_. 
This is similar to a cycling-type measurement and is 
quantum non-demolition (after wait time k). There will 
be limitations impose on size of £ due to finite anhar- 
monicity of the TCQ. 

To achieve tuning of g+ we modify the original cir- 
cuit and replace the Josephson junctions by SQUIDs 
with Josephson energy E^ and E^ (this is hinted 
at in Fig. 1 B). In making this replacement the only 
change in the above theory is the replacement Ej ± — » 

£™ ax cos(7r^ ± /^>o)v /l + ^ 2 tan 2 (7T<^/$ ) with £™ ax = 
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FIG. 3. (color online) Matrix element of the collective Cooper 
pair number operator (A) and transition energy (B) of the 
dark state as a function of the energy ratio Ej + /Ec for 
n g+ = n' g _ — 0, Ei — —Ec and Ej_ is numerically solved to 
ensure that only the coupling strength is tuned (blue) and fre- 
quency (red) . Solid lines are from a numerical diagonalization 
and dashed lines are from the coupled anharmonic oscillator 
model. 



E j± + e j± > d = ( e j1 ~ e j1 ) l E J± ax and $ *± is the exter " 
nal flux applied to each SQUID which we assume to be 

independent (this is not required but simplifies our argu- 
ment). This independent control allows us to change Ej ± 
independently which in-turn allows independent control 
on g + and u q . To illustrate this we consider the sym- 
metric case and plot in Fig. 3 the normalized coupling 
strength g+h/2e 2 V rms f3 (A) and uj q (B) as a function of 
the ratio Ej + /Ec when Ej = — Ec and Ej_ is numeri- 
cally solved to ensure that only the coupling rate (blue) 
and frequency (red) vary respectively for both the full nu- 
merical (solid) and effective model (dashed). In the full 
numerical model g+ = 2e 2 V rms + /3_n_) |0)/ft. 

Here the independent control is clearly observed. Note 
that while our numerical investigation was only for the 
symmetric case independent tunable g+ (from zero to 
large values) and uo q will still occur when the device is 
not symmetric. There is just a different condition on 
Ej ± for the required tuning. 

With the extra control channel there is the possibility 
of additional qubit decoherence from flux fluctuations. 
A reasonable estimate for the flux induced relaxation is 
7i ux = E± \(0\dH/d$ x± \l)\ 2 MlS I± (w q )/h 2 where M± 
is the mutual inductance between the bias line and the 
TCQ, and Si ± (uj q ) is the current noise in the bias line 
which at low temperatures Si ± (uj q ) ~ hu q /R [27]. Tak- 
ing M = 200 $ /A R = 50 ft, EJ^/h = 20 GHz, 
E c /h = Ejjh = 0.35 GHz, and d = 10% we find 7\ w Is. 
To estimate the contribution to dephasing we assume 
the noise is 1/f and from Ref. [28] the dephasing time 
« \duj q /d$> x± \~ 1 /A ( f ) where is the flux noise mea- 
sured at 1 Hz which for similar superconducting devices 
has been measured to be 10 _6 ^>o / V^Hz [29] . This predicts 
a « 20/is (at Q x± = 7r$o/4). This rates are consis- 
tent with the T\ and predictions for the transmon 



(with only frequency control) from flux noise [7]. That 
is, in comparison to the transmon the additional flux line 
that gives both independent control of the coupling and 
frequency has not added any extra noise. 

In conclusion we have presented a new device for quan- 
tum information processing with superconducting cir- 
cuits. It is an extension of the transmon and uses quan- 
tum interference to achieve independent control of both 
the coupling strength and frequency. Furthermore it can 
be tuned to a configuration where it is in a decoherence 
free subspace with respect to relaxation induced by the 
Purcell effect whilst still allowing efficient readout of its 
quantum state. It also offers the possibility of imple- 
menting a cycling- type measurement. 
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